Accumulation of Microswimmers due to Their Collisions with a Surface 
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In this letter we propose a kinematic model to show how collisions with a surface and rotational 
Brownian motion give rise to the accumulation of micro-swimmers near a surface. In this model, an 
elongated microswimmer invariably travels parallel to the surface after hitting it from any incident 
angle. It then swims away from the surface after some time, facilitated by rotational Brownian 
motion. Simulations based on this model reproduce the density distributions measured for the 
small bacteria E. coli and Caulohacter crescentus, as well as for the much larger bull spermatozoa 
swimming in confinement. 



Swimming aids the function and development of mi- 
croorganisms in many ways. For example, it enhances 
the formation of biofilms, which provide favorable mi- 
croenvironments for bacteria to cope with environmental 
stresses jH- Swimming also helps transport spermato- 
zoa toward eggs for fertilization Interestingly, the 
number density of cells as a function of distance from a 
surface has been measured for E. coli Q and bull sperma- 
tozoa [3|, showing in both cases values much higher near 
the surface than far away. This near surface accumula- 
tion has mainly been attributed to a hydrodynamic at- 
traction between the cells and the surface 0,01. Recently, 
Berke et al. Q combined the effects of the hydrodynamic 
attraction and the translational Brownian motion of the 
cells to predict the distribution of E. coli as a function of 
distance. As noted by the authors however, this in- 
terpretation is not applicable to cells within 10 /xm from 
the surface, where most accumulation occurs. 

In this letter we present a different account for the near 
surface accumulation. We ignore the hydrodynamic at- 
traction but emphasize the role of the collision with a 
surface at low Reynolds number and rotational Brown- 
ian motion in a confined environment. We show that a 
typical microswimmer with an elongated shape will swim 
parallel to a surface after hitting it and therefore accu- 
mulate near the surface. Rotational Brownian motion Q 
then relaxes the accumulation by randomly changing the 
swimming direction so that the cells have chances to swim 
away from the surface. In the extreme case of no rota- 
tional Brownian motion, all the cells would end up swim- 
ming in close proximity with the surface. In the oppo- 
site extreme of very fast rotational Brownian motion, the 
cells will quickly change to any possible swimming di- 
rection and subsequently would be found anywhere with 
equal probability. In reality, a microorganism randomly 
changes its swimming direction with a finite rotational 
diffusion constant, resulting in a distribution in between 
the two extremes, that is, more cells stay near the surface 
and fewer far away. 

We used the bacterium C. crescentus strain CB15 
SB3860, which is a CheR mutant of Apilin (YB375) 
and swims forward exclusively, to examine the details 
of near surface swimming. Swarmer cells of this mu- 
tant do not follow circular trajectories when swimming 
near surfaces The strains were synchronized with 



the plate releasing method 0, 0] to obtain cultures with 
primarily swimming cells. The synchronized culture was 
sealed between a glass slide and a coverslip with vacuum 
grease for optical microscopy observation. Broken cov- 
erslip pieces were used as spacers so that the thickness 
of the microscopy sample is ^--^200 /im. A 20x objective 
(Nikon Plan Apo, NA 0.75) was used on a Nikon E800 
microscope to take 5 snap shots of swimming cells at 0.1 
second intervals under darkfield mode using a CoolSnap 
CCD camera (Princeton Instruments) and MetaMorph 
software (Universal Imaging). The objective was focused 
on planes 5, 15, 25 jim, etc., away from the surface and 
the cell number distribution was measured following the 
method of Berke et al. Q . We noted that although this 
objective has a 1.4 ^m depth of field, cells up to nearly 5 
/im off the focal plane appeared as bright spots. There- 
fore the measured cell density was an average over a ~10 
fim thick layer. 

The swimming speed and the rotational diffusion con- 
stant were obtained from the videos taken for cells over 
50 fim away from both surfaces. The average swimming 
speed was ^45 ^m/s. The rotational diffusion constant 
was measured from ~200 swimming trajectories. The 
swimming direction at moment t was taken as the di- 
rection from the position at t to the position at t + 0.1 
s. With this definition the change in direction Aip over 
time interval At was obtained and the rotational diffu- 
sion constant Dr was calculated to be 0.12 rad^/s, using 
the equation < Aip'^ >= 2DrAt. 

With particular interest we examined 3-D trajectories 
as the cells approached and swam near a surface, until 
they took off. To do so, we focused the objective on 
the top surface and recorded the swimming trajectories. 
Example trajectories are shown in Fig. [1^ by overlaying 
consecutive frames taken at the rate of 10 frames per 
second. The cell body appeared as a sharp bright spot 
when it was in the focal plane and as a ring when it 
was away. Wu et al. [§| found that the ring size was 
proportional to the distance of the cell from the focal 
plane and therefore can be calibrated to determine the 
distance. Two examples of 3-D trajectories of the cells 
1 and 2 in Fig. [1^ are plotted in Figs. [TJd and Ic. Most 
cells approached the surface at an angle and then swam 
parallel to the surface for some time before leaving. The 
manner of C. crescentus hitting a surface is similar to 
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Figure 1: Trajectories of Caulobacter swarmer cells swimming 
near a glass surface, (a) Overlay of consecutive darkfield im- 
ages taken at 10 frames per second, (b) and (c) are 3-D plots 
(red) and projections (blue) on the glass surface of the trajec- 
tories of cells 1 and 2 in (a). Arrows indicate the swimming 
directions. 
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Figure 2: (a) Force and torque analysis of a forward swimming 
cell hitting a surface, (b) Required time (solid) for the cell to 
become parallel to the surface and the rotation rate (dashed) 
as function of angle 6. 



that of E. coli observed with three dimensional tracking 
microscopy [lo| . 

We analyzed the force and torque on C. crescentus 
swimming near a surface and found that it would invari- 
abley swim parallel to the surface shortly after hitting 
the surface. For purposes of the model, we approximate 
the cell as a sphere attached with a helical filament of 
length L propelled by a longitudinal force Fp. After the 
cell hits the surface at an angle 9, its velocity component 
along the direction normal to the surface (y-axis Fig.[2t) 
becomes zero. It will have a swimming speed Vx along 
the X-axis and a rotation rate along the z-axis (not 
shown in the figure) . We ignore the increase in hydrody- 
namic drag on the cell due to the nearby surface 0, [lH 
and assume that the surface only provides a force Fs to 
stop the swimming along the y-axis. The hydrodynamic 
drag forces on the whole cell (sphere plus helical filament) 
are split into components parallel and perpendicular to 
the long axis, Fy and F±. The hydrodynamic torque T 
on the whole cell is depicted with respect to the sphere 
center. The forces and torque are given by 
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where V]| = cos 9 and V± = Vx sin 9 are the speed 
components along and perpendicular to the helical axis 
and A is the friction matrix, where Aij > and A23 = 
A32. 

At a low Reynolds number, the force balance along x- 
axis is Fp cos 9+F^^ cos 9+F± sin = and torque balance 
along 2-direction is F = 0, which together determine the 



swimming speed and rotation rate as 
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Since A22A33 > A23, the denominator in the expressions 
above are always positive. In the case as shown in Fig. [2^, 
Vx > and > 0. Therefore the cell swims toward the 
right and the filament rotates toward the surface. 

We can estimate how fast the cell turns parallel to 
the surface during the collision. Mathematically, the cell 
would take an infinitely long time to become parallel to 
the surface, as calculated from Eq. [31 In practice, how- 
ever, since the rotational Brownian motion of C. crescen- 
tus varies its orientation by 0.1 rad within less than 0.1 
sec, we estimate instead the time needed for the cell ahgn- 
ment with the surface to fall under 0.1 rad. The parame- 
ters for a typical C. crescentus are An = 2.2x10^® 
Nsm-i, A22 = 2.5 X 10"8 Nsm'S A33 = 1.9 x lO^^^ 
Nms, and A23 ~ 5.3 x 10~^^ Ns. The propulsive force is 
Fp = AiiV ^ 1 X 10-12 N, where V is the bulk swim- 
ming speed. The rotation rate after hitting a surface is 
shown in Fig. [IJa, which reaches 9 rad/s at 55°. If a cell 
hits the surface at an angle 9o, the time for it to become 
parallel is d9/n (Fig. [lb). This is less than 0.2 s for a 
typical angle of 6*0 = 30°, and less than 0.3 s for an angle 
as large as 85°. Therefore in the following discussion we 
state in a practical sense that a cell becomes parallel to 
the surface after the collision. 
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Figure 3: (a) Rod model of a microswimmer near surface. 
The black end represents that of the cell body and the gray 
end the flagellar filament, (b) An example of simulated dis- 
tance h (red) and angle (p (blue) as functions of time for the 
microswimmer, using the parameters of C. crescentus. The 
two surfaces are separated by 200 /im. 



Now we examine how a swimming microorganism takes 
off after hitting a surface. To further simpUfy the model, 
we approximate the elongated swimmer propelled by a 
longitudinal force as a nonuniform rod (Fig. [SK). This 
rod swims forward at speed V in the bulk fluid. The rod 
has a rotation center at position O, which is of a distance 
Li away from the head and L2 away from the tail. The 
head has a larger drag per unit length than the tail does, 
and thus Li < L2. Due to the small size, the rod under- 
goes constant Brownian motion with a rotational diffu- 
sion constant Dr and translational diffusion constant Dt. 
Since An ^ A22, we ignore the angle dependence of Dt. 

The change in distance of the rotation center to the 
surface y is determined by the translational Brown- 
ian motion and the swimming direction, which is con- 
stantly altered by the rotational Brownian motion. Over 
a time interval At, Ay — X^sin^At + Q^2DtAt, and 
A(j) = <^V2£'rAt, where ^ and <^ are random numbers 
with zero mean and unit variance. The translational 
Brownian motion contributes much less than swimming 
to the displacement for microorganisms swimming at tens 
of ^m/s. When near the surface, the changes in distance 
and angle are also restricted by the solid surface to sat- 
isfy y ^ Li sin(— 0) when the head is closer to the surface 
and y ^ L2 sin when the tail is closer. Similar restric- 
tions hold when a cell is near the top surface. Knowing 
Dt and Dr, we can track the distance y and angle (j) over 
time. The distance of the head from the surface h, which 
is what was measured in the experiments, is determined 
hj h = y + Li sin (j). The probability distribution of a cell 
at distance h is obtained by tracking a cell swimming 
between the two surfaces over 10^ - 10^ sec. 



We simulated the distance and angle of swimming C. 
crescentus between two glass surfaces separated by 200 
fim. The cell was treated as a i = 6 fim rod, with a typ- 
ical translational diffusion constant on the order of 0.1 
/im^/s and the measured rotational diffusion constant of 
0.12 rad^/s. The rotation center was approximated at a 
position where Li = 0.3L. Fig.jSb shows examples of dis- 
tance (red) and angle (blue) varying over time. The cell 
hits the top and bottom surfaces repeatedly as it swims 
between them. The simulated distance from the bottom 
surface was recorded every 0.1 second and a histogram of 
distances was made using a bin size of 10 /im. The simu- 
lated distribution is plotted in Fig. [4] (blue) and compared 
with the measured one for C. crescentus (triangle). The 
simulation clearly shows higher densities near the sur- 
faces, with the entire proflle in excellent agreement with 
the measurements. 

This model is also applicable to the distribution of 
E. coli and bull spermatozoa between two surfaces. We 
took the cell number distribution of E. coli from refer- 
ence and that of bull spermatozoa from reference 01, 
converted them to probability density, and plotted them 
in comparison with that of C. crescentus in Fig.jH E. coli 
is similar in size to C. crescentus and it is not surprising 
that they have similar distributions. Bull spermatozoa 
are ten times larger yet surprisingly the distribution is 
similar to that of bacteria. Nevertheless, this similarity 
is actually predicted by our model. To simulate for bull 
spermatozoa, we treated it as a 60 ^m long rod with 
a translational diffusion constant on the order of 0.01 
/im^/s and a rotational diffusion constant of Dr = 10^'' 
rad^/s, which is ~1000 times smaller than that of C. 
crescentus. The simulation shows that the difference in 
distribution between the bull spermatozoa (red) and the 
C. crescentus (red) is so small that it cannot be distin- 
guished by the observations under the set conditions. 

The distribution of swimming cells near a surface can 
be understood intuitively based on the variation of angle 
due to rotational Brownian motion. Once a cell hits a 
surface, the angle becomes zero and the cell swims par- 
allel to the surface for some time. The cell then swims 
away from the surface by changing swimming direction 
due to rotational Brownian motion. Let us examine, for 
instance, the time ti and ^2 needed for a cell to swim to 
different distances from the surface, hi < h2. Since the 
cell swims to distance hi before reaching ^-2,^1 < ^2 must 
hold. Assuming the swimming direction is only affected 
by rotational Brownian motion, the mean square angle 
at time t is < cjP >~ 2Drt. Therefore, the closer a cell 
swims near the surface, the smaller the angle is. The per- 
pendicular component of swimming speed Vy = V sin 
is smaller when closer to the surface. The dwell time of 
a cell staying at a distance h is inversely proportional 
to Vy and therefore the probability density is larger near 
the surface. Additional features of the distribution can 
be understood based on this physical picture. When the 
cell swims close to the surface after a coUision, </) <C 1. 
The probability density is ~ 1/0, which drops sharply 
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Figure 4: Comparison between simulated density distribu- 
tions at rotational diffusion constants 0.12 (blue) and 0.0001 
(red) rad^/s and the measured distributions of C. crescentus 
(up triangle), E. coli (down triangle, ref [3|), and bull sperma- 
tozoa (square, ref [3|). The separation between the two glass 
surfaces is 200 ^tm. Inset compares simulated distribution at 
rotational diffusion constants of 10 (diamond), 1 (circle), 0.1 
(triangle), and 0.0001 (square) rad^/s at a swimming speed of 
50 fi,m/s, corresponding to rod lengths of ~1.3,~2.8, ~6, and 
~60 nm, respectively. The dotted lines indicate the probabil- 
ity density if there is no surface accumulation. 

with the increasing angle and distance. This picture also 
predicts that a microswimmer with a large rotational 
diffusion constant leaves a surface more rapidly after a 
collision and hence the accumulation near the surface is 
weaker, as shown for Dr = 10 rad^/s in the inset of Fig.lH 

In this model we have ignored the hydrodynamic in- 
teraction between the cell and the surface. In reality, 
the swimming cell generates a flow, which interacts with 
the nearby surface, reorienting and attracting the cell to- 
wards the surface. Berke et al. [3] calculated this effect 



for bacteria when the cell is >10 /xm away from the sur- 
face. A simple estimation shows that this effect is small 
compared to rotational Brownian motion and swimming 
when the cell is nearly parallel to the surface. For ex- 
ample, at a distance h — \0 /im and angle (j) = 0-1 rad, 
calculation based on their model yields a reorientation 
rate of '--^O.Ol rad/s and an attraction speed of ~1 /^m/s, 
while in 1 second the rotational Brownian motion can re- 
orient the cell by 0.4 rad on average and the component 
of swimming speed normal to the surface is on the order 
of 10 ^m/s. Therefore, the modification to the distri- 
bution due to the long range hydrodynamic interaction 
is small when the distance is >10 ^m. When a cell is 
less than 1 ^m from the surface, the large hydrodynamic 
friction between the cell and surface may keep the cells 
near the surface for a long time [llj. The effect of this 
extension of dwell time is not that dramatic since the 
distribution is binned by 10 ^m in distance. The hydro- 
dynamic interaction in the range of 1 to 10 ^m is yet to 
be described theoretically. Its effect on distribution of 
cells in this range, however, is expected to be secondary 
as evident by the good agreement between the simulation 
results ignoring it and the experimental measurements. 

In conclusion, we have demonstrated the effects of col- 
lision and rotational Brownian motion on swimming mi- 
croorganisms near surfaces. The collision with a surface 
resets the swimming direction to be parallel to it and 
the rotational Brownian motion then randomly alters 
the swimming direction, which leads to the accumula- 
tion. An excellent agreement is obtained between the 
simulations based on this picture and the experimental 
results. Recently, various artificial microswimmers have 
been developed |13l . [H, , which usually have an elon- 
gated shape. Similar effect of the collision and rotational 
Brownian motion is expected when they swim near a sur- 
face. 
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